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Abstract 



In the spectrum of systems showing chaos-assisted tunneling, three-state 
crossings are formed when a chaotic singlet intersects a tunnel doublet. We 
study the dissipative quantum dynamics in the vicinity of such crossings. A 
harmonically driven double well coupled to a bath serves as a model. Markov 
and rotating-wave approximations are introduced with respect to the Floquet 
spectrum of the time-dependent central system. The resulting master equa- 
tion is integrated numerically. We find various types of transient tunneling, 
determined by the relation of the level width to the inherent energy scales of 
the crossing. The decay of coherent tunneling can be significantly retarded 
or accelerated. Modifications of the quantum asymptotic state by the cross- 
ing are also studied. The comparison with a simple three-state model shows 
that in contrast to the undamped case, the participation of states outside the 
crossing cannot be neglected in the presence of dissipation. 
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I. INTRODUCTION 



In a quantum setting, the coexistence of regular and chaotic regions in a mixed phase 
space leads to a variety of uncommon coherence phenomena. Most prominent among them 
is chaotic tunneling the coherent exchange of probability between symmetry-related 

regular islands that are separated by a chaotic layer, not by a static potential barrier. Chaotic 
tunneling comes about by an interplay of classical nonlinear — typically bistable — dynamics 
and quantum coherence. It therefore reflects features of the classical phase space, such as the 
width of the chaotic layer []7||| and the structure of the "coast line" Jll],|l2j separating it from 



the adjacent regular regions, as well as fine details of the quantum spectrum like exact and 
avoided crossings. In previous work |7|J§, it has been shown that, as the chaotic layer grows 
with increasing nonlinearity, the tunnel splittings widen collectively. Superimposed on this 
global trend of the parameter dependence, however, there are strong fluctuations, occurring 
on a smaller parameter scale and restricted to individual tunnel doublets |J, IT, IS, 13]. A 
major source of these fluctuations are the disturbances of doublets, suffered as they are 
intersected by other levels. The most common type of such intersections is formed when 
a doublet, pertaining to a pair of eigenstates located on symmetry- related quantizing tori, 



encounters a singlet that belongs to an eigenstate in the chaotic sea [0. The chaotic states, 



even if they no longer come in pairs close in energy, can still be classified as even or odd. This 
fact determines the structure of the singlet-doublet crossings: The partner in the doublet 
sharing the symmetry of the chaotic singlet is repelled by the singlet — together, they form 
an avoided crossing. The partner with opposite symmetry must either be intersected in 
an exact crossing close to the avoided one, or else the order within the doublet is reversed 
from one side of the crossing to the other. The variety of configurations of such crossings is 
sketched in Fig. |l[ 

In the present paper, we study chaotic tunneling in the vicinity of singlet-doublet cross- 
ings, under the influence of incoherent processes. Near a crossing, level separations deviate 
vastly, in both directions, from the typical tunnel splitting (cf. Fig. [I]). This is reflected 
in time-domain phenomena ranging from the suppression of tunneling to a strong increase 
in its rate and to complicated quantum beats |TIJ. In Section |Q[ we briefly review chaotic 
tunneling in our working model, a harmonically driven quartic double well |15|-[T7|. Singlet- 
doublet crossings are identified and characterized by their signature in terms of quasienergy 
and mean energy. We describe the behaviour of the eigenstates close to a crossing and 
analyze the coherent dynamics in terms of a simple three-state model. 

Tunneling is associated with extremely small energy scales, all the more in the semiclas- 
sical regime we are interested in. It is therefore particularly sensitive to any disruption of 
coherence as it occurs due to the unavoidable coupling to the environment. As immediate 
consequences, the symmetry underlying the formation of tunnel doublets is generally bro- 
ken, and an additional energy scale is introduced, the effective finite width attained by each 
discrete level. Tunneling and related coherence phenomena are thus rendered transients that 
occur — if at all — on the way towards an asymptotic equilibrium state. Singlet-doublet cross- 
ings, in turn, drastically change the nondissipative energy scales and replace the two-level 
by a three- level structure. As a consequence, the familiar way tunneling fades out in the 
presence of dissipation is also significantly altered. Near a crossing, the coherent dynamics 
can last much longer than for the unperturbed doublet, despite the presence of the same 
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decoherence than outside the crossing, establishing "chaos-induced coherence" . Depending 
on temperature, it can also be destroyed on a much shorter time scale. 

In Section [TIT], we first outline the modifications of our model necessary to incorporate 
dissipation. We treat the quantum dissipative dynamics within the framework of the Floquet 



formalism [|18| |23|1 . In the present context of an interplay of tunneling with incoherent 
processes, we focus on the regime of weak dissipation. Very long time scales are therefore 
involved, both in the coherent dynamics and in the bath response. In this case, the usual 
(Markov and rotating-wave) approximations resorted to in the elimination of the external 
degrees of freedom, have to be reconsidered critically. The coupling to the environment 
indirectly couples the three levels in the crossing to all the other states of the central system. 
On the basis of numerical results for the full driven double well with dissipation, we reveal 
the limitations of the three-level approximation and identify additional features of the full 
dynamics not covered by it. In particular, we consider the long-time asymptotics — the 
quantum attractor — and the phase-space structure associated with it. Section |V| serves to 
summarize our results and to suggest directions of further research. 



II. CONSERVATIVE CLASSICAL AND QUANTUM DYNAMICS 

A. The model and its symmetries 

As a prototypical working model, we consider the quartic double well with a spatially 
homogeneous driving force harmonic in time. It is defined by the Hamiltonian 

Hxm{x,p;t) = H Q (x,p) + H^t), 

ffo(*, P )=^4* 2 + ^, a) 

Hi(x;t) = Sxcos(u>t). 

Apart from mere scaling, the classical phase space of H (x, p) depends only on the presence or 
absence, and the sign, of the x 2 term. Otherwise, it has no free parameter. In the quantum- 
mechanical case, the parameter D/h plays the role of an inverse quantum of action. It 
controls the barrier height and can be interpreted as the (approximate) number of doublets 
with energies below the top of the barrier. Accordingly, the classical limit is reached by 
letting D/h — > oo. The influence of the driving on the classical phase space structure is 
characterized by the rescaled amplitude 

F = S/V8D (2) 

and frequency to. This implies that the classical dynamics is independent of the barrier 
height D. 

The unperturbed Hamiltonian Ho(x,p) is invariant under the parity P: x — > —x, p — > —p, 
t — > t. This symmetry is generally destroyed by the driving. With the above choice of 



Hi(x;t), however, a more general, dynamical symmetry is retained |15|-[l7|,^J . It is defined 
by the operation 

P w : x — > —x, p — > —p, t^t + n/uj, (3) 
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and represents a generalized parity acting in the extended phase space spanned by x, p, and 
phase, i.e., time tmod (2ir/u). The invariance under of H(x,p; t) allows to classify all its 
quantum eigenstates as even or odd. 

As a consequence of the periodic time dependence of H(x,p;t), the relevant generator 
of the quantum dynamics is now the Floquet operator p5|-[29|1 

U = Texp(--j o dtH(t)\, (4) 

where T denotes time ordering. According to the Floquet theorem, the adiabatic states of 
the system are the eigenstates of U. They can be written in the form 

\Mt))=e- i€at/h \<t>a(t)), (5) 

with 



\(j> a {t + 2ir/u)) = \(/> a (t)). 

Expanded in the basis spanned by these Floquet states, the propagator of the driven system 
reads 

U(t',t)=Y / e- i ^ t '- t)/h \Ut'))(Mt)\- (6) 

a 

The associated eigenphases e a , referred to as quasienergies, come in classes, e a , n = e a + nhuj, 
n = 0, ±1, ±2, . . .. This is suggested by a Fourier expansion of the \4> a (t)), 



-iniot 



\4>oc{t)) = \ C »,n)e 
n 

(jj r2n/w 

Kn) = 7T dt\<f) a (t))e m ^. (7) 

Z7T JO 

The index n counts the number of quanta in the driving field. Otherwise the members of a 
class a are physically equivalent. Therefore, the quasienergy spectrum can be reduced to a 
single "Brillouin zone", —hu/2 < e < huj/2. 

Since the quasienergies have the character of phases, they can be ordered only locally, 
not globally. A quantity that is defined on the full real axis and therefore does allow for a 



complete ordering, is the mean energy [p8|,p£ 

Ea = %- r /UJ dt(i; a (t)\H(t) \i> a (t)) = ((Mt)\H(t) |0 Q (t)». (8) 

Z7T JO 

It is related to the corresponding quasienergy by 

E a = e a + (((f> a (t)\m^\(l>a(t))), (9) 

where the outer angle brackets denote the time average over one period of the driving, as 
indicated in eq. (pi). The second term on the right-hand side plays the role of a geometric 
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phase accumulated over this period pOfl . Without the driving, E a = e a , as it should be. By 
inserting the Fourier expansion (|7|) the mean energy takes the form 



E a = J2(e a + nhu) {c a>n \c at n) . 



(10) 



It shows that the n th Floquet channel gives a contribution e a + nhu to the mean energy, 
weighted by the Fourier coefficient {c a ^ n \c a , n ) . 

Quasienergies and Floquet states are obtained numerically by solving the matrix eigen- 
value equation p5|,[28| 



n' k> 



(11) 



equivalent to the time-dependent Schrodinger equation. It is derived by inserting the eigen- 
states (|5|) into the Schrodinger equation, Fourier expanding, and using the representation 
by the eigenstates of the unperturbed Hamiltonian, H Q \^k) = Ek\^k)- We introduced the 
abbreviations 

H n ,k;n',k' = (E k — nflLj)5 n ^ n i5k-k' 

+ - S Xk t k' {$n-l-ri + $n+l~n>), 

Cn,k (^fel^n)) 
Xk,k> = (^k\x |*fc/). 

The invariance of the system under is of considerable help in solving Eq. (jTT|) , because 
it completely decouples the respective systems of eigenvalue equations for the even and odd 
subspaces 



B. Coherent chaotic tunneling 

With the driving Hi(x;t) switched off, the classical phase space generated by H(x,p;t) 
exhibits the constituent features of a bistable Hamiltonian system. There is a separatrix at 
E = 0. It forms the border between two sets of trajectories: One set, with E < 0, comes in 
symmetry-related pairs, each partner of which oscillates in either one of the two potential 
minima. The other set consists of unpaired trajectories, with E > 0, that encircle both wells 
in a spatially symmetric fashion. 

Increasing the amplitude of the driving from zero onwards has two principal consequences 
for the classical dynamics (Fig. |^): The separatrix is destroyed as a closed curve and replaced 
by a homoclinic tangle |3l| of stable and unstable manifolds. As a whole, it forms a chaotic 



layer in the vicinity and with the topology of the former separatrix. It opens the way for 
diffusive transport between the two potential wells. Due to the nonlinearity of the potential, 
there is an infinite set of resonances of the driving with the unperturbed motion, both 
inside and outside the wells Since the period of the unperturbed, closed trajectories 



diverges for E — > 0, the resonances accumulate towards the separatrix of the unperturbed 
system. By its sheer phase-space area, the first resonance (the one for which the periods of 
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the driving and of the unperturbed oscillation are in a ratio of 1:1) is prominent among the 
others and soon (in terms of increasing amplitude F) exceeds the size of the "order-zero" 
regular areas near the bottom of each well. For the values of F and uo chosen in the numerical 
calculations in this paper, all higher resonances remain negligible in size. The borderline 
between the chaotic layer along the former separatrix and the regular regions within and 
outside the wells is therefore quite sharply defined. The "coastal strip" formed by hierarchies 
of regular islands around higher resonances remains narrow. For the tunneling dynamics, 
the role of states located in the border region |TT is therefore not significant here. 

Both major tendencies in the evolution of the classical phase space — extension of the 
chaotic layer and growth of the first resonance — leave their specific traces in the quasienergy 
spectrum. The tunnel doublets characterizing the unperturbed spectrum for E < pertain 
to states located on pairs of symmetry-related quantizing tori in the regular regions within 
the wells. With increasing size of the chaotic layer, the quantizing tori successively resolve in 
the chaotic sea. The corresponding doublets disappear as distinct structures in the spectrum 
as they attain a splitting of the same order as the mean level separation. The gradual 
widening of the doublets proceeds as a smooth function of the driving amplitude |7||§. 



This function roughly obeys a power law [34], cf. Fig. H. As soon as a pair of states is no 



longer supported by any torus-like manifold, including fractal [35 and vague tori |55fl , the 
corresponding eigenvalues detach themselves from the regular ladder to which they formerly 
belonged. They can then fluctuate freely in the spectrum and thereby "collide" with other 
chaotic singlets or regular doublets. 

The appearance of a regular region, large enough to accomodate several eigenstates, 
around the first resonance introduces a second ladder of doublets into the spectrum. Size 
and shape of the first resonance vary in a way different from the fate of the main regular 
region. The corresponding doublet ladder therefore moves in the spectrum independently of 
the doublets that pertain to the main regular region, and of the chaotic singlets. This gives 
rise to additional singlet-doublet and even to doublet-doublet encounters. 



C. Three-level crossings 

Among the various types of quasienergy crossings that occur according to the above 
scenario, those involving a regular doublet and a chaotic singlet are the most common. In 
order to give a quantitative account of such crossings and the associated coherent dynamics, 
and for later reference in the context of the incoherent dynamics, we shall now discuss them 
in terms of a simple three-state model, devised much in the spirit of Ref . || . 

Far to the left of the crossing, we expect the following situation: There is a doublet of 
Floquet states 

l^ + (t))=e-^|0 r + (t)), (12) 
|^(t))=e-^ ++A ^|0 r -(t)), (13) 

with even (superscript '+') and odd ('— ') generalized parity, respectively, residing on a pair 
of quantizing tori in one of the regular (subscript 'r') regions. We have assumed that the 
quasienergy splitting (as opposed to the unperturbed splitting) is — e+ = A > 0. The 
global relative phases can be chosen such that the superpositions 
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liM*)> = 73 (lVv + (*)> ± IVC(*)>) (14) 

are initially localized in the right and the left well, respectively, and tunnel back and forth 
with a frequency A/% given by the tunnel splitting in the presence of the driving. 
As the third player, we introduce a Floquet state 

lC(t)> = e- lfe++A+Ao)t/R |C(t)>, (15) 

located mainly in the chaotic (subscript 'c') layer, so that its time-periodic part |$T(t)) 
contains a large number of harmonics. Without loss of generality, its generalized parity is 
fixed to be odd. For the quasienergy, we have assumed that e~ = e+ + A + A c , where |A C | 
can be regarded as a measure of the distance from the crossing. 

The structure of the classical phase space then implies that the mean energy of the 
chaotic state should be close to the top of the barrier and far above that of the doublet. We 
assume, like for the quasienergies, a small splitting of the mean energies pertaining to the 
regular doublet, E~ — <C E~ — Ef. 

In order to model an avoided crossing between \ip~) and |^), we suppose that there is 
a non- vanishing matrix element ((<f>~\Hjyw\<f)~)) = b > 0. For the singlet- doublet crossings 
under study, we typically find that A <C b <C Tiuj. Neglecting the coupling with all other 
states, we obtain the three-state (subscript '3s') Floquet Hamiltonian 

/0 \ 
Hss = 4 + A b , (16) 
\0 b A + A c / 

in the three-dimensional Hilbert space spanned by {\<f>^(t)), \4>~(t)), \<f)~(t))}. Its Floquet 
states read 

|Vo + (*)> = e-^i^rc*)), 

l^r(*)> = ^ t/% (\<P7(t)) cos/3 - \<p-(t)) sin/3) , (17) 
IV£ (*)> = ^ t/h (|0T(*)> sm/? + |CW> cos/3) . 
Their quasienergies are 

4 = 4, el 2 = et + A+ l -A cT \^Al + Ab^ (18) 
The mean energies are approximately given by 

^o" = E 4 

E^ = E~ cos 2 + E~ sin 2 /3, (19) 
E 2 = E~ sin 2 (3 + E~ cos 2 /3, 

where contributions of the matrix element b have been neglected. The angle (3 describes the 
mixing between the Floquet states and and is a measure of the distance to the 
avoided crossing. By diagonalizing the Hamiltonian (|TB|), we obtain 
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2(3 = arctan 



'26 
A, 



< f3 < 



71 



(20) 



For (3 — > 7r/2, corresponding to — A c ^> 6, we retain the situation far left of the crossing, 
as outlined above, with ~ \^c)i IV^) ~ To the far right of the crossing, i.e., for 

(3 — > or A c ^> 6, the exact eigenstates ) and |^) have interchanged their identity with 
respect to the phase-space structure |]ID| . Here, we have ) ~ |^~) and ~ |V\T)- The 
mean energy is essentially determined by the phase-space structure. Therefore, there is also 
an exchange of E± and E 2 in an exact crossing, cf. Eq. (|TP|), while .Eg" remains unaffected 
(Fig. ^b). The quasienergies and must intersect close to the avoided crossing of e{ 
and e 2 (Fig. |]a). Their crossing is exact, since they pertain to states with opposite parity 
(cf. Fig.|a,b). 

Numerical evidence shows that this idealized picture is not always correct. It may well 
happen that even far away from a crossing, the doublet splitting does not exactly return to 
its value on the opposite side (cf. Figs. |5]a, ||). It is even possible that an exact crossing of 
€q and €1 does not take place at all in the vicinity of the crossing. In that case, the relation 
of the quasienergies in the doublet gets reversed via the crossing (Fig. |l]c,d). Nevertheless, 
the above scenario captures the essential features. 

To study the dynamics of the tunneling process, we focus on the state 



\m) 



V2 V 



A4t/% \<t>t{t)) +e- i ^|0rW)cos/3 + e- i ^^|0 2 -(t))sin/5) . 



(21) 



It is constructed such that at t — 0, it corresponds to the decomposition of (\ip*) + \ip~)) /y/2 
(cf. Eq. ([14])) in the basis (|17D at finite distance from the crossing. Therefore, it is initially 
localized in the regular region in the right well and follows the time evolution under the 
Hamiltonian fll6[). From Eqs. fll4|), (|17D, we find the probabilities for its evolving into |^r), 
IV'l), or j^c), respectively, to be 



p K (t) = \(Mt)\m)f 



1 / Eit 2 e 2 t . o „ 

- 1 + cos — — cos (3 + cos — sin p + 

2 \ n n 



p L (t) = \(Mt)\m)f 



2 a e ^t . 2 a , 

cos — cos p — cos — sin p + 
n n 



cos ■ 



cos • 



(ei - e 2 )t 



(fi - e 2 )t 

h 



cos 2 (3 sin 2 (3 



cos 2 (3 sin 2 (3 



p c (t) = immm 



cos • 



(ei - e 2 )t 
h 



cos 2 (3 sin 2 f3. 



(22) 



In sufficient distance from the crossing, there is only little mixing between the regular 
and the chaotic states, i.e., sin/3 <C 1 or cos/5 <C 1. The tunneling process then follows 
the familiar two-state dynamics involving only an d iV'r"); with tunnel frequency A/h 
(Fig. 0a). 

Close to the avoided crossing, cos (3 and sin (3 are of the same order of magnitude, and 
\i[>x), \ip2~) become very similar to one another. Both have now support in the chaotic layer 
as well as in the symmetry-related regular regions and thus are of a hybrid nature. Here, the 
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tunneling involves all the three states and must at least be described by a three-level system. 
The exchange of probability between the two regular regions proceeds via a "stop-over" in 
the chaotic region [H,[K]]. The three quasienergy differences that determine the time scales 
of this process are in general all different, leading to complicated beats (Fig. |7|b). 

However, for A c = —A/2, the two quasienergies ej~ and — degenerate. At this point, 
which marks the center of the crossing, the number of different frequencies in the three-level 
dynamics reduces to two again. This restores the familiar coherent tunneling in the sense 
that there is again a simple periodic exchange of probability between the regular regions 
IIJ . However, the rate is much larger if compared to the situation far off the crossing, and 



the chaotic region is now temporarily populated during each probability transfer, twice per 
tunneling cycle (Fig. 0c). 

In order to illustrate the above three-state model and to demonstrate its adequacy, we 
have numerically studied a singlet-doublet crossing that occurs for the double-well potential, 
Eq. (0), with D/fi = 4, at a driving frequency uo = 0.982 and amplitude F = S/V8D = 
0.015029 (Figs. [|, ||). A comparison (not shown) of the appropriately scaled three-state 
theory (Fig. |j) with this real singlet- doublet crossing (Fig. |5]) shows satisfactory agreement. 
Note that in the real crossing, the quasienergy of the chaotic singlet decreases as a function 
of F, so that the exact crossing occurs to the left of the avoided one. 



III. INCOHERENT QUANTUM DYNAMICS 
A. Master equation 

1. System-bath model 

To achieve a microscopic model of dissipation, we couple the system (|TJ) bilinearly to a 
bath of non-interacting harmonic oscillators |J7| . The total Hamiltonian of system and bath 
is then given by 



oo 



"M^W + E^ + ^-^j j. (23) 

We here couple the position x of the system to an ensemble of oscillators with masses m u , 
frequencies u u , momenta p u , and coordinates x u , with the coupling strength g u . The bath 
is fully characterized by the spectral density of the coupling energy, 

Jn = nf-^-5(uj-^). (24) 

For the time evolution we choose an initial condition of the Feynman- Vernon type: at t — t , 
the bath is in thermal equilibrium and uncorrelated to the system, i.e., 

P{to) = Ps(*o) ® PB,eq, (25) 

where pe.eq = exp(— PH B )/tr B exp(—(3H B ) is the canonical ensemble of the bath and 1//3 = 
k B T. 
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Due to the bilinearity of the system-bath coupling, one can always eliminate the bath 
variables to get an exact, closed integro-differential equation for the reduced density matrix 
Ps(t) = tr B p(t), which describes the dynamics of the central system, subject to dissipation 



2. Born-Markov approximation 

In most cases, however, the integro-differential equation for ps(t) can be solved only 
approximately. In particular, in the limit of weak coupling, 

7 < k B T/h, (26) 
7<k-ea'IA (27) 

it is possible to truncate the time-dependent perturbation expansion in the system-bath 
interaction after the second-order term. The quantity 7, to be defined below, denotes the 
effective damping of the dissipative system, and \e a — e a >\/% are the transition frequencies 
of the central system. In the present case, the central system is understood to include 
the driving p0|-f2"3|j , so that the transition frequencies are given by quasienergy differences. 



The autocorrelations of the bath decay on a time scale h/k^T and thus in the present limit, 
instantaneously on the time scale I/7 of the system correlations. With the initial preparation 
(p5|) , the equation of motion for the reduced density matrix in this approximation is given 
by@ 



i 1 f°° 

p s {t) = -- [H s (t),p s (t)} + — / du J(u)n th {hu) 

tl 71 tl J —00 

x dr (e iujT [x(t - r, t)p s {t),x] + H.c.) (28) 

where x(t',t) denotes the position operator in the interaction picture defined by 

x{t',t) = U\t',t)xU{t\t), (29) 

with U(t',t), the propagator of the conservative driven double well, given in Eq. (||). 'H.c' 
denotes the Hermitian conjugate and 

nth ^ = e e/k B l _ I = -^th(-e) - 1 (30) 

is the thermal occupation of the bath oscillator with energy e. To achieve a more compact 
notation, we require J(— to) = —J(oj). In the following, we shall restrict ourselves to an 
Ohmic bath, J(u>) = m'juj. This defines the effective damping constant 7. 

We use the time-periodic components \4> a {t)) of the Floquet states as a basis to expand 
the density operator, Eq. (p8l). Expressing the matrix elements 

Xvpit) = (</> a (t)\x\Mt)) ( 31 ) 

of the position operator by their Fourier coefficients 
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X a p, n = ((U^xe-^lMt))) = X; a ,- n , (32) 
X a p(t) = e muJt X a p^ n , (33) 

n 

yields for the matrix elements a a p of the reduced density matrix ps the equation of motion 



M i' L'|.l>:; 



cr a/ 3(t) = f t (Mt)\Ps(t)\Mt)) 

+ X {Naa',nX aa ', n O'a.'P'Xpip n i 
ce'P'nn' 

— N a ' /3' , n X aa ' , n 'X aa ' )n a/3' ^e 1 ^ 11 ^ + H.c. . (34) 

Note that the coefficients of this differential equation are periodic in time with the period 
of the driving. The N a p >n are given by 

N a( 3, n = N(e a - 6p + nfku), N(e) = —^-n^e). (35) 

For e 3> k^T, N(e) approaches zero. 

Since the position operator x is odd under transformation with the generalized parity @, 
the master equations (^) and (|34| ) are invariant under transformation with P w . Therefore, 
both the conservative and the dissipative dynamics preserve the parity of the operator 
1 0a) (0/3 1- It is even if |0 Q ) and \4>p) belong to the same parity class and odd otherwise. 
Note that in particular, the projectors |0 a )(0 a | and thus all density matrices diagonal in the 
Floquet basis are even under P w . 



3. Rotating-wave approximation 
Assuming that dissipative effects are relevant only on a time scale much larger than the 



period 2tc/uj of the driving, we average the coefficients of the master equation fl34|) over 2-k/uj 
and end up with the equation of motion 



&ap{t) — — *( e « — e p)Vap{t) + X/ L a (3a> fi'Va' f3> , (36) 
n a'P' 

with the time-independent dissipative part 

L a pa'P' = X {Xaa',n + Nppt >n ) X aa / jn Xpip^ n 



n 



~$PP' X/ ^ 'P" ' a ' ,nX a p" - n Xpn a i >n 
P",n 

—0~aa' X N a "P',nXp> a »-nX a »p tn . (37) 
a"n 

This step amounts effectively to a rotating-wave approximation. It is, however, less restric- 
tive than the rotating-wave approximation introduced in [p0|,21]1 where dissipative effects are 



averaged over the generally longer time scale m&x aj p :Tl (2irh/(e a — ep + nhco)). 
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B. Dissipative chaos-assisted tunneling 



The crucial effect of dissipation on a quantum system is the disruption of coherence: 
a coherent superposition evolves into an incoherent mixture. Thus, phenomena based on 
coherence, such as tunneling, are rendered transients that fade out on a finite time scale 
^decoh- In general, for driven tunneling in the weakly damped regime, this time scale gets 



shorter for higher temperatures, reflecting the growth of transition rates ||39|| . However, 
there exist counterintuitive effects. For example, in the vicinity of an exact crossing of the 
ground-state quasienergies, the coherent suppression of tunneling |l5j can be stabilized with 



higher temperatures |T5|,ff5LBD| until levels outside the doublet start to play a role. We have 



studied dissipative chaos-assisted tunneling numerically, using again the real singlet-doublet 
crossing introduced in Sect. |11 C] (cf. Fig. ^|) as our working example. 

In the vicinity of a singlet- doublet crossing, the tunnel splitting increases significantly — 
the essence of chaos-assisted tunneling. During the tunneling, the chaotic singlet becomes 
populated periodically with frequency |ej — ejf cf. Eq. ( |22|) and Fig. [?]. The high mean 
energy of this singlet results in an enhanced decay of coherence at times when \i/j c ) is popu- 
lated (Fig. |8|). For the relaxation towards the asymptotic state, also the slower transitions 
within doublets are relevant. Therefore, the corresponding time scale t re i ax can be much 
larger than t docoh (Fig. g). 

To obtain quantitative estimates for the dissipative time scales, we approximate tdecoh 
by the decay rate of trp 2 , a measure of coherence, averaged over time, 

-Jo"**™ < 38) 
i(trp 2 (0)-trp 2 ( f )). (39) 



^decoh 



Because of the stepwise decay of the coherence (Fig. |8|), we have chosen the propagation time 
t as an n-fold multiple of the duration 27r/|e^ — | of the chaotic beats. For this procedure 
to be meaningful, n should be so large that the coherence decays substantially during the 
time t (in our numerical studies to a value of approximatly 0.9). The time scale t re iax of the 
approach to the asymptotic state is given by the reciprocal of the smallest real part of the 
eigenvalues of the dissipative kernel. 

Outside the singlet-doublet crossing we find that the decay of coherence and the relax- 
ation take place on roughly the same time scale (Fig. [Hj). At F « 0.013, the chaot ic singlet 
induces an exact crossing of the ground-state quasienergies (see Fig. ^J), resulting in a sta- 
bilization of coherence with increasing temperature. At the center of the avoided crossing, 
the decay of coherence becomes much faster and is essentially independent of temperature. 
This indicates that transitions from states with mean energy far above the ground state play 
a crucial role. 



C. Asymptotic state 

As the dynamics described by the master equation (|28|) is dissipative, it converges in 
the long-time limit to an asymptotic state Poo(i). In general, this attractor remains time 
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dependent but shares all the symmetries of the central system, i.e. here, periodicity and 
generalized parity. However, the coefficients of the master equation (|37|) for the matrix 
elements <j a p, valid within a rotating- wave approximation, are time independent and so the 
asymptotic solution also is. This means that we have eliminated the explicit time dependence 
of the attractor by representing it in the Floquet basis and introducing a mild rotating-wave 
approximation. 

To gain some qualitative insight into the asymptotic solution, we focus on the diagonal 
elements 

Laaa'a' = 2 N aa > <n \ X aa i <n \ 2 , a ^ a', (40) 

n 

of the dissipative kernel. They give the rates of the direct transitions from \<p a ') to \(f> a ). 



Within a cruder rotating-wave approximation ||21|| , these are the only non- vanishing con- 
tributions to the master equation which affect the diagonal elements a aa of the density 
matrix. 

In the case of zero driving amplitude, the Floquet states \<p a ) reduce to the eigenstates 
of the undriven Hamiltonian i^DW- The only non- vanishing Fourier component is then 
\c a} o), and the quasienergies e a reduce to the corresponding eigenenergies E a . Thus L aaa i a i 
consists of a single term proportional to N(e a — e a >) only. It describes two kinds of thermal 
transitions: decay to states with lower energy and, if the energy difference is less than k^T, 
thermal activation to states with higher energy. The ratio of the direct transitions forth and 
back then reads 

Laaa'a' I \^a ^a')\ i , ai\ 

1 = exp V^F~ ■ (41) 

J^a'a'aa \ ^B 1 / 

We have detailed balance. Therefore, the steady-state solution o aa ' (oo) ~ 
exp(— Ea/ksT) 5 a a'- In particular, the occupation probability decays monotonically with 
the energy of the eigenstates. In the limit k^T —>■ 0, the system tends to occupy the ground 
state only. 

For a strong driving, each Floquet state \<f) a ) contains a large number of Fourier compo- 
nents and Laaa'a' is given by a sum over contributions with quasienergies e a — e a ' + nhuj. 
Thus a decay to states with "higher" quasienergy (recall that quasienergies do not allow 
for a global ordering) becomes possible due to terms with n < 0. Physically, they describe 
dissipative transitions under absorption of driving-field quanta. Correspondingly, the sys- 
tem tends to occupy Floquet states comprising many Fourier components with low index n. 
According to Eq. (|TUD, these states have low mean energy. 



The effects under study are found for a driving with a frequency of the order of unity. 
Thus for a quasienergy doublet, i.e., far off the three-level crossing, we have |e Q — e a > \ hu, 
and La'a'aa is dominated by contributions with n < 0, where the splitting has only small 
influence. However, as a consequence of symmetry, the splitting is the main difference 
between the two partners of the quasienergy doublet. Therefore, with respect to dissipation, 
both should behave similarly. In particular, one expects an equal population of the doublets 
even in the limit of zero temperature (Fig. |TT|a). This is in contrast to the undriven case. 

In the vicinity of a singlet- doublet crossing the situation is more subtle. Here, the odd 
partner, say, of the doublet mixes with a chaotic singlet, cf. Eq. (|TT|), and thus acquires 
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components with higher energy. Due to the high mean energy E~ of the chaotic singlet, 
close to the top of the barrier, the decay back to the ground state can also proceed indirectly 
via other states with mean energy below E~ . Thus \<pi) and \<f>2 ) are depleted and mainly 
|0o~) will be populated. However, if the temperature is significantly above the splitting of the 
avoided crossing, thermal activation from \4>q) to |0^~ 2 ), accompanied by depletion via the 
states below E~ , becomes possible. Thus asymptotically, all these states become populated 
in a steady flow (Figs. [H]b,c). 

An important global characteristic of the asymptotic state is its coherence trp^,. Its 
value gives approximately the reciprocal of the number of incoherently occupied states. It 
equals unity only if the attractor is a pure state. According to the above scenario, we 
expect tr to assume the value 1/2, in a regime with strong driving but preserved doublet 
structure, reflecting the incoherent population of the ground-state doublet. In the vicinity 
of the singlet-doublet crossing where the doublet structure is dissolved, its value should be 
close to unity for temperatures k^T <C b and much less than unity for k-^T ^> b (Figs. O, 



13|) . This means that the crossing of the chaotic singlet with the regular doublet leads to an 



improvement of coherence if the temperature is below the splitting of the avoided crossing, 
and a loss of coherence for temperatures above the splitting. This phenomenon amounts to 
a chaos-induced coherence or incoherence, respectively. 

The crucial role of the decay via states not involved in the three-level crossing can 
be demonstrated by comparing it with the dissipative dynamics including only these three 
levels (plus the bath). At the crossing, the three-state model results in a completely different 
type of asymptotic state (Fig. |l^). The failure of the three-state model in the presence of 
dissipation clearly indicates that in the vicinity of the singlet- doublet crossing, it is important 
to take a large set of levels into account. 



IV. CONCLUSION 

Nonlinear systems with a mixed phase space still represent a formidable challenge for 
a theoretical understanding. On the classical level, the intricate fractal interweaving of 
regular and chaotic regions of phase space is far from being exhaustively studied. Quantum- 
mechanical uncertainty reduces the richness of the classical phase-space structure, but at the 
same time, coherence effects like tunneling add new elements to the dynamics. For example, 
in a bistable system, a driving with an energy in the regime between the unperturbed tunnel 
splittings and the typical separation of doublets leads to transport phenomena that combine 
tunneling and classical chaotic diffusion as essential elements. Yet another basic energy scale 
enters if the unavoidable coupling of the system to its environment is taken into account. 
If the finite width the levels attain in the presence of dissipation is comparable, in turn, to 
the energies characterizing driven tunneling, unfamiliar phenomena due to the interplay of 
chaos, tunneling, and decoherence occur in the transient dynamics and possibly also in the 
asymptotic state of the driven damped quantum system. 

In this paper, we have selected a specific case out of the host of problems to be studied 
in this field: dissipative chaotic tunneling in the vicinity of crossings of chaotic singlets with 
tunnel doublets. In order to obtain a firm quantitative basis, we have endowed a prototypical 
model for chaotic tunneling, a harmonically driven double well, with dissipation. Thereby, 
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we have followed the usual approach of coupling to a heat bath, but adapted to the periodic 
time dependence of the central system. Moreover, we have carefully avoided to destroy, by 
the approximations introduced in the derivation of a master equation, the specific spectral 
characteristics of chaotic tunneling that we are interested in. As a simple intuitive model to 
compare against, we have constructed a three-state system which in the case of vanishing 
dissipation, provides a faithful description of an isolated singlet-doublet crossing. 

As an example for the interplay of chaotic tunneling with dissipation, we mention an 
effect that might be termed "chaos- induced coherence" for short. Suppose that damping 
and temperature are such that far off the singlet-doublet crossing, the level width is of the 
same order or larger than the tunnel splitting. There, coherent tunneling is then largely 
suppressed even as a transient. Close to the crossing, however, the widening of the doublet 
enforced by the intersecting third level creates a separation of the time scales of the coherent 
tunneling and of its decay due to incoherent processes, without altering the dissipation as 
such. As a consequence, the transient dynamics now exhibits many tunneling cycles, each 
of which includes two "stop overs" in the chaotic phase-space region, one on the way to, one 
on the way back. 

The study of the asymptotic state, the "quantum attractor" , demonstrates clearly that 
a three-state model of the singlet-doublet crossing is no longer adequate once dissipation is 
effective. This is so because the coupling to the heat bath enables processes of decay and 
thermal activation that connect the states in the crossing with other, "external" states of the 
central system. In the presence of driving, the asymptotic state is no longer literally a state 
of equilibrium. Rather, incoherent processes create a steady flow of probability involving 
states within as well as outside the crossing. As a result, the composition of the asymptotic 
state, expressed for example by its coherence trp^,, will be markedly different at the center 
of the crossing as compared to the asymptotic state far away from the crossing, even if that 
is barely visible in the corresponding phase-space structure. 

Many more phenomena at the overlap of chaos, tunneling, and dissipation await being 
unraveled. They include four-state crossings formed when two doublets intersect, chaotic 
Bloch tunneling along extended potentials with a large number of unit cells instead of just 



two, the influence of decoherence on a multi-step mechanism of chaotic tunneling [TT|JTj] and 
transient tunneling between coexisting strange attractors, to mention only a few. At the 
same time, analytical tools more specifically taylored to the investigation of the quantum- 



classical correspondence in mixed systems, such as a phase-space entropy will provide 
additional insight. 
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FIGURES 




FIG. 1. Possible configurations of quasienergy crossings between a chaotic singlet and a regular 



doublet. Different line types signify different parity. See Sect. [II C| for the labeling of the levels. 
Note that only for configurations a,b, the order of the regular doublet is restored in passing through 
the crossing. In configurations c,d, it is reversed. 
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FIG. 2. Stroboscopic classical phase-space portraits, at uit = 2irn, of the harmonically driven 
quartic double well, Eq. (|), at F = 0.015, uj = 0.982. 
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FIG. 3. Dependence on the driving amplitude of the splitting of the ground-state doublet, of 
the harmonically driven quartic double well, Eq. ([!]), at D/fi = 8 and u> = 0.95. Diamonds are 
quasienergy, crosses mean-energy splittings. The dips interrupting the smooth parameter depen- 
dence indicate crossings of the ground-state doublet with chaotic singlets. After ||. 
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FIG. 4. A singlet-doublet crossing, according to a three-state model, Eq. ([if]), in terms of the 
dependence of the quasienergies (panel a) and the mean energies (b) on the coupling parameter 
A c /b. Unperturbed energies are marked by dotted lines, the energies for the case with coupling by 
full lines for even and dashed lines for odd states. 
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FIG. 5. Singlet-doublet crossing found numerically for the driven double well, Eq. (|l|), at 
D/h = 4 and oj = 0.982, in terms of the dependence of the quasienergies (a) and the mean energies 
(b) on the driving amplitude F = S/\/8D. Values of the driving amplitude used in Figs. are 
marked by dotted vertical lines. Full and dashed lines indicate energies of even and odd states, 
respectively. Bold lines give the mean energies of the chaotic singlet and the ground-state doublet 
depicted in panel a. 
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FIG. 6. Splitting of the lowest doublets for D/h = 4 and cj = 0.982. The arrows indicate the 
locations of the exact and the avoided crossing within a three-level crossing of the type sketched 
in Fig. [l]a. 




FIG. 7. Time evolution of a state initially localized in the right well, in the vicinity of the 
singlet-doublet crossing shown in Fig. [|, in terms of the probabilities to be in the right well (return 
probability, marked by full lines), in the reflected state in the left well (dashed), or in the chaotic 
state |V> C > (dotted). Parameter values are as in Fig. g, and F = 0.0145 (a), 0.0149 (b), 0.015029 
(c). 
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FIG. 8. Occupation probabilities as in Fig. |7|a,c, but in the presence of dissipation. The 
dash-dotted line shows the time evolution of tr p 2 . The parameter values are D/h = 4, u = 0.982, 
ry = 10" 6 , k B T/h = 10" 4 , and F = 0.0145 (a), 0.015029 (b). The inset in (a) is a blow up of the 
rectangle in the upper left corner of that panel. 
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FIG. 9. Time evolution of the return probability Pr (full line) and the coherence function 
tr p 2 (dash-dotted) during loss and regain of coherence. The parameter values are as in Fig. ||]b. 
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FIG. 10. Time scales of the decay of tr p 2 (a) and of the relaxation towards the asymptotic 
solution (b) near the singlet-doublet crossing. Near the exact crossing (F ~ 0.013, full vertical 
line) coherence is stabilized, whereas at the center of the avoided crossing (F « 0.015, dashed 
vertical line) the decay of coherence is accelerated. The parameter values are D/h = 4, u = 0.982, 
7 = 10 -6 , temperature as given in the legend. 
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FIG. 11. Occupation probability a aa of the Floquet states \ip a ) in the long-time limit. The 
parameter values are D/h = 4, u = 0.982, 7 = 10" 6 , and F = 0.013 (a), 0.0145 (b), 0.015029 (c), 
temperature as given in the legend. 
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FIG. 12. Coherence of the asymptotic state in the vicinity of a singlet-doublet crossing for 
different temperatures as given in the legend. The other parameter values are D/K = 4, to = 0.982, 
and 7 = 1(T 6 . 




FIG. 13. Coherence of the asymptotic state in the vicinity of a singlet-doublet crossing for 
F = 0.013 (a) and F = 0.015029 (b): exact calculation (full line) compared to the values resulting 
from a three-level description (dashed) of the dissipative dynamics. The other parameter values 
are D/h = 4, u = 0.982, and 7 = 10~ 6 . 
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